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Abstract 

This article provides an introduction to Schramm(stochastic)-Loewner evolution 
(SLE) and to its connection with conformal field theory, from the point of view 
of its application to two-dimensional critical behaviour. The emphasis is on the 
conceptual ideas rather than rigorous proofs. 
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1 Introduction 

1.1 Historical overview 

The study of critical phenomena has long been a breeding ground for new ideas in the- 
oretical physics. Such behaviour is characterised by a diverging correlation length and 
cannot easily be approximated by considering small systems with only a few degrees of 
freedom. Initially, it appeared that the analytic study of such problems was a hopeless 
task, although self-consistent approaches such as mean field theory were often successful 
in providing a semi-quantitative description. 

Following Onsager's calculation of the free energy of the square lattice Ising model in 
1944, steady progress was made in the exact solution of an ever-increasing number of 
lattice models in two dimensions. While many of these arc physically relevant, and the 
techniques used have spawned numerous spin-offs in the theory of integrable systems, it 
is fair to say that these methods have not cast much light on the general nature of the 
critical state. In addition, virtually no progress has been made from this direction in 
finding analytic solutions to such simple and important problems as percolation. 

An important breakthrough occurred in the late 1960's, with the development of renor- 
malisation group (RG) ideas by Wilson and others. The fundamental realisation was 
that, in the scaling limit where both the correlation length and all other macroscopic 
length scales are much larger than that of the microscopic interactions, classical critical 
systems are equivalent to renormalisable quantum field theories in euclidean space-time. 
Since there is often only a finite or a dcnumerablc set of such field theories with given 
symmetries, at a stroke this explained the observed phenomenon of universality: systems 
with very different constituents and microscopic interactions nevertheless exhibit the same 
critical behaviour in the scaling limit. This single idea has led to a remarkable unification 
of the theoretical bases of particle physics, statistical mechanics and condensed matter 
theory, and has led to extensive cross-fertilisation between these disciplines. These days, 
a typical paper using the ideas and methods of quantum field theory is as likely to appear 
in a condensed matter physics journal as in a particle physics publication (although there 
seems to be a considerable degree of conservatism among the writers of field theory text 
books in recognising this fact.) 

Two important examples of this interdisciplinary flow were the development of lattice 
gauge theories in particle physics, and the application of conformal field theory (CFT), 
first developed as a tool in string theory, to statistical mechanics and condensed matter 
physics. As will be explained later, in two-dimensional classical systems and quantum 
systems in 1-1-1 dimensions conformal symmetry is extremely powerful, and has led to a 
cornucopia of new exact results. Essentially, the RG programme of classifying all suitable 
renormalisable quantum field theories in two dimensions has been carried through to its 
conclusion in many cases, providing exact expressions for critical exponents, correlation 
functions, and other universal quantities. However the geometrical, as opposed to the 
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algebraic, aspects of conformal symmetry are not apparent in this approach. 

One minor but nevertheless theoretically influential prediction of these methods was the 
conjectured crossing formula for the probability that, in critical percolation, a cluster 
should exist which spans between two disjoint segments of the boundary of some simply 
connected region (a more detailed account of this problem will be given later.) With this 
result, the simmering unease that mathematicians felt about these methods came to the 
surface (see, for example, the comments in 0.) What exactly are these renormalised local 
operators whose correlation functions the fleld theorists so happily manipulate, according 
to rules that sometimes seem to be a matter of cultural convention rather than any 
rigorous logic? What does conformal symmetry really mean? Exactly which object is 
conformally invariant? And so on. Aside from these deep concerns, there was perhaps 
also the territorial feeling that percolation theory, in particular, is a branch of probability 
theory, and should be understood from that point of view, not merely as a by-product of 
quantum fleld theory. 

Thus it was that a number of pure mathematicians, versed in the methods of probability 
theory, stochastic analysis and conformal mapping theory, attacked this problem. Instead 
of trying make rigorous the notions of fleld theory about local operators, they focused 
on the random curves which form the boundaries of clusters on the lattice, and on what 
should be the properties of the measure on such curves in the continuum limit as the lattice 
spacing approaches zero. The idea of thinking about lattice models this way was not 
new: in particular in the 1980s it led to the very successful but non-rigorous Coulomb gas 
approach [B' to two-dimensional critical behaviour, whose results parallel and complement 
those of CFT. However, the new approach focused on the properties of a single such curve, 
conditioned to start at the boundary of the domain, in the background of all the others. 
This leads to a very speciflc and physically clear notion of conformal invariance. Moreover, 
it was shown by LoewnerJ^ in the 1920s that any such curve in the plane which does not 
cross itself can be described by a dynamical process called Loewner evolution, in which the 
curve is imagined to be grown in a continuous fashion. Instead of describing this process 
directly, Loewner considered the evolution of the analytic function which conformally 
maps the region outside the curve into a standard domain. This evolution, and therefore 
the curve itself, turns out to be completely determined by a real continuous function 
at- For random curves, at itself is random. (The notation at is used rather than a{t) 
to conform to standard usage in the case when it is a stochastic variable.) Schramm j3] 
argued that, if the measure on the curve is to be conformally invariant in the precise sense 
referred to above, the only possibility is that at be one-dimensional Brownian motion, with 
only a single parameter left undetermined, namely the diffusion constant k. This leads to 
stochastic-, or Schramm-, Loewner evolution (SLE). (In the original papers by Schramm 
et al. the term 'stochastic' was used. However in the subsequent literature the 'S' has 
often been taken to stand for Schramm in recognition of his contribution.) It should 
apply to any critical statistical mechanics model in which it is possible to identify these 
non-crossing paths on the lattice, as long as their continuum limits obey the underlying 
conformal invariance property. For only a few cases, including percolation, has it been 
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proved that this property holds, but it is beheved to be true for suitably defined curves 
in a whole class of systems known as 0(n) models. Special cases, apart from percolation, 
include the Ising model. Potts models, the XY model, and self-avoiding walks. They each 
correspond to a particular choice of k. 

Starting from the assumption that SLE describes such a single curve in one of these 
systems, many properties, such as the values of many of the critical exponents, as well as 
the crossing formula mentioned above, have been rigorously derived in a brilliant series 
of papers by Lawler, Schramm and Werner (LSW)|ni. Together with Smirnov's proof [7j 
of the conformal invariance property for the continuum limit of site percolation on the 
triangular lattice, they give a rigorous derivation [S] of the values of the critical exponents 
for two-dimensional percolation. This represents a paradigm shift in rigorous statistical 
mechanics, in that results are now being derived directly in the continuum for models for 
which the traditional lattice methods have, so far, failed. 

However, from the point of view of theoretical physics, these advances are perhaps not so 
important for being rigorous, as for the new light they throw on the nature of the critical 
state, and on conformal field theory. In the CFT of the 0(n) model, the point where a 
random curve hits the boundary corresponds to the insertion of a local operator which 
has a particularly simple property: its correlation functions satisfy linear second-order 
differential equations^. These equations turn out to be directly related to the Fokker- 
Planck type equations one gets from the Brownian process which drives SLE. Thus there 
is a close connection, at least at an operational level, between CFT and SLE. This has 
been made explicit in a series of papers by Bauer and Bernard ^U] (see also fl].) Other 
fundamental concepts of CFT, such as the central charge c, have their equivalence in 
SLE. This is a rapidly advancing subject, and some of the more recent directions will be 
mentioned in the concluding section of this article. 

1.2 Aims of this article 

The original papers on SLE are mostly both long and difficult, using, moreover, concepts 
and methods foreign to most theoretical physicists. There are reviews, in particular those 
by Werner and by Lawler |13j which cover much of the important material in the 
original papers. These are however written for mathematicians. A more recent review by 
Kager and NienhuisjMJ describes some of the mathematics in those papers in way more 
accessible to theoretical physicists, and should be essential reading for any reader who 
wants then to tackle the mathematical literature. A complete bibliography up to 2003 
appears in fT^. 

However, the aims of the present article are more modest. First, it does not claim to 
be a thorough review, but rather a semi-pedagogical introduction. In fact some of the 
material, presenting some of the existing results from a slightly different, and hopefully 
clearer, point of view, has not appeared before in print. The article is directed at the 
theoretical physicist familiar with the basic concepts of quantum field theory and critical 
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behaviour at the level of a standard graduate textbook, and with a theoretical physicist's 
knowledge of conformal mappings and stochastic processes. It is not the purpose to prove 
anything, but rather to describe the concepts and methods of SLE, to relate them to other 
ideas in theoretical physics, in particular CFT, and to illustrate them with a few simple 
computations, which, however, will be presented in a thoroughly non-rigorous manner. 
Thus, this review is most definitely not for mathematicians interested in learning about 
SLE, who will no doubt cringe at the lack of preciseness in some of the arguments and 
perhaps be puzzled by the particular choice of material. The notation used will be that of 
theoretical physics, for example (■ ■ ■) for expectation value, and so will the terminology. 
The word 'martingale' has just made its only appearance. Perhaps the largest omission 
is any account of the central arguments of LSWjHI which relate SLE to various aspects of 
Brownian motion and thus allow for the direct computation of many critical exponents. 
These methods are in fact related to two-dimensional quantum gravity, whose role in this 
is already the subject of a recent long article by Duplantier[16j. 
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2 Random curves and lattice models. 
2.1 The Ising and percolation models 

In this section, we introduce the lattice models which can be interpreted in terms of 
random non-intersecting paths on the lattice whose continuum limit will be described by 
SLE. 

The prototype is the Ising model. It is most easily realised on a honeycomb lattice (see 
Fig. P). At each site r is an Ising 'spin' s(r) which takes the values ±1. The partition 




Figure 1: Ising model on the honeycomb lattice, with loops corresponding to a term in the 
expansion of (^. Alternatively, these may be thought of as domain walls of an Ising model on 
the dual triangular lattice. 

function is 

Zising = Tr e^^S..'^««(-') oc Tr n (1 + xs{r)s{r')) , (1) 

rr' 

where x = tanh f3J, and the sum and product are over all edges joining nearest neighbour 
pairs of sites. The trace operation is defined as Tr = Hr (| J2s{r) )y so that Trs(r)" = 1 if 
n is even, and if it is odd. 

At high temperatures [jSJ <^ 1) the spins are disordered, and their correlations decay 
exponentially fast, while at low temperatures (/5J ^ 1) there is long-range order: if 
the spins on the boundary are fixed say, to the value then (s(r)) ^ even in the 
infinite volume limit. In between, there is a critical point. The conventional approach 
to the Ising model focuses on the behaviour of the correlation functions of the spins. In 
the scaling limit, they become local operators in a quantum field theory (QFT). Their 
correlations are power-law behaved at the critical point, which corresponds to a massless 
QFT, that is a conformal field theory (CFT). From this point of view (as well as exact 
lattice calculations) it is found that correlation functions like (s(ri)s(r2)) decay at large 
separations according to power laws |ri — r2|~^^: one of the aims of the theory is to obtain 
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the values of the exponents x as well as to compute, for example, correlators depending 
on more than two points. 

However, there is an alternative way of thinking about the partition function (^, as 
follows: imagine expanding out the product to obtain 2^ terms, where N is the total 
number of edges. Each term may be represented by a subset of edges, or graph Q, on 
the lattice, in which, if the term xs{r)s{r') is chosen, the corresponding edge (rr') is 
included in Q, otherwise it is not. Each site r has either 0, 1, 2 or 3 edges in Q. The 
trace over s(r) gives 1 if this number is even, and if it is odd. Each surviving graph 
is then the union of non- intersecting closed loops (see Fig. Q). In addition, there can be 
open paths beginning and ending at a boundary. For the time being, we suppress these 
by imposing 'free' boundary conditions, summing over the spins on the boundary. The 
partition function is then 

■Rising = ^ X , (2) 

g 

where the length is the total of all the loops in Q. When x is small, the mean length of 
a single loop is small. The critical point Xc is signalled by a divergence of this quantity. 
The low-temperature phase corresponds to x > Xc- While in this phase the Ising spins 
are ordered, and their connected correlation functions decay exponentially, the loop gas is 
in fact still critical, in that, for example, the probability that two points lie on the same 
loop has a power-law dependence on their separation. This is the dense phase. 

The loops in Q may be viewed in another way: as domain walls for another Ising model 
on the dual lattice, which is a triangular lattice whose sites R lie at the centres of the 
hexagons of the honeycomb lattice (see Fig. Q). If the corresponding interaction strength 
of this dual Ising model is (/?</)*, then the Boltzmann weight for creating a segment of 
domain wall is e"^*-'^"'^*. This should be equated to x = tanh(/5J) above. Thus we see that 
the high-temperature regime of the dual model corresponds to low temperature in the 
original model, and vice versa. Infinite temperature in the dual model {{jSJ)* = 0) means 
that the dual Ising spins are independent random variables. If we colour each dual site 
with s{R) = +1 black, and white if s{R) = —1, we have the problem of site percolation 
on the triangular lattice, critical because Pc = ^ for that problem. Thus the curves with 
X = 1 correspond to percolation cluster boundaries. (In fact in the scaling limit this is 
believed to be true throughout the dense phase x > Xc-) 

So far we have discussed only closed loops. Consider the spin-spin correlation function 

/ / X / N\ Tr s{n)s{r2) Urr' (1 + xs{r)s{r')) 
= Trn..(l + a:.(r).(rO) ' 

where the sites ri and r2 lie on the boundary. Expanding out as before, we see that the 
surviving graphs in the numerator each have a single edge coming into ri and r2. There 
is therefore a single open path 7 connecting these points on the boundary (which does 
not intersect itself nor any of the closed loops.) In terms of the dual variables, such a 
single open curve may be realised by specifying the spins s{R) on all the dual sites on 
the boundary to be +1 on the part of the boundary between ri and r2 (going clockwise) 
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Figure 2: The exploration process for tlie Ising model. At each step the walk turns L or R 
according to the value of the spin in front of it. The relative probabilities are determined by the 
expectation value of this spin given the fixed spins either side of the walk up to this time. The 
walk never crosses itself and never gets trapped. 

and —1 on the remainder. There is then a single domain wall connecting ri to r2. SLE 
describes the continuum limit of such a curve 7. 

Note that we could also choose r2 to lie in the interior. The continuum limit of such 
curves is then described by radial SLE fSec. 13. 

2.1.1 Exploration process 

An important property of the ensemble of curves 7 on the lattice is that, instead of 
generating a configuration of all the s{R) and then identifying the curve, it may be 
constructed step-by-step as follows (see Fig. [21). Starting from ri, at the next step it should 
turn R or L according to whether the spin in front of it is +1 or —1. For independent 
percolation, the probability of either event is |, but for x < 1 it depends on the values of 
the spins on the boundary. Proceeding like this, the curve will grow, with all the dual sites 
on its immediate left taking the value +1, and those its right the value —1. The relative 
probabilities of the path turning R or L at a given step depend on the expectation value 
of the spin on the site R immediately in front of it, given the values of the spins already 
determined, that is, given the path up to that point. Thus the relative probabilities that 
the path turns R or L are completely determined by the domain and the path up to that 
point. This implies the crucial 

Property 2.1 (lattice version). Let 71 be the part of the total path covered after a 
certain number of steps. Then the conditional probability distribution of the remaining 
part of the curve, given 71, is the same as the unconditional distribution of a whole curve, 
starting at the tip and ending at r2 in the domain P \ 71. 

In the Ising model, for example, if we already know part of the domain wall, the rest 
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of it can be considered as a complete domain wall in a new region in which the left and 
right sides of the existing part form part of the boundary. This means the path is a 
history-dependent random walk. It can be seen (Fig. |2|) that when the growing tip r 
approaches an earlier section of the path, it must always turn away from it: the tip never 
gets trapped. There is always at least one path on the lattice from the tip r to the final 
point r2. 



2.2 0(n) model 

The loop gas picture of the Ising and percolation models may simply be generalised by 
counting each closed loop with a fugacity n: 

rv \ ^ length number of loops 

^O(n) — 2^ X n 

s 

This is called the 0(n) model, for the reason that it gives the partition function for 
n-component spins s(r) = (si(r), . . . Sn{r)) with 

^o(n) = Trn(l + ^s(r)-s(r')) , (5) 

rr' 

where Tr Sa{r)sb{r) = Sab- Following the same procedure as before we obtain the same set 
of closed loops (and open paths) except that, on summing over the last spin in each closed 
loop, we get a factor n. The model is called 0(n) because of its symmetry under rotations 
of the spins. The version makes sense only when n is a positive integer (and note that 
the form of the partition function is different from that of the conventional 0(?t,) model, 
where the second term is exponentiated.) The form in Q is valid for general values of 
n, and it gives a probability measure on the loop gas for real n > 0. However, the dual 
picture is useful only for n = 1 and n = 2 (see below.) As for the case n = 1, there is a 
critical value Xc{n) at which the mean loop length diverges. Beyond this, there is a dense 
phase. 

Apart from n = 1, other important physical values of n are: 



(4) 



• n = 2. In this case we can view each loop as being oriented in either a clockwise 
or anti-clockwise sense, giving it an overall weight 2. Each loop configuration then 
corresponds to a configuration of integer valued height variables h{R) on the dual 
lattice, with the convention that the nearest neighbour difference h{R') — h{R) takes 
the values 0, +1 or —1 according to whether the edge crossed by RR' is unoccupied, 
occupied by an edge oriented at 90° to RR', or at —90°. (That is, the current running 
around each loop is the lattice curl of h.) The variables h{R) may be pictured as 
the local height of a crystal surface. In the low-temperature phase (small x) the 
surface is smooth: fluctuations of the height differences decay exponentially with 
separation. In the high-temperature phase it is rough: they grow logarithmically. In 
between is a roughening transition. It is believed that relaxing the above restriction 
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on the height difference h{R') — h{R) does not change the universahty class, as long 
as large values of this difference are suppressed, for example using the weighting 
exp [ — i3{h{R') — h{R))'^]. This is the discrete Gaussian model. It is dual to a model 
of 2-component spins with 0(2) symmetry called the XY model. 

• n = 0. In this case, closed loops are completely suppressed, and we have a single 
non- self-intersecting path connecting ri and r2, weighted by its length. Thus, all 
paths of the same length are counted with equal weight. This is the self-avoiding 
walk problem, which is supposed to describe the behaviour of long flexible polymer 
chains. As x — > Xc—, the mean length diverges. The region x > Xc is the dense 
phase, corresponding to a long polymer whose length is of the order of the area of 
the box, so that it has finite density. 

• n = —2 corresponds to the loop-erased random walk. This is an ordinary random 
walk in which every loop, once it is formed, is erased. Taking n — —2 in the 0{n) 
model of non-intersecting loops has this effect. 

2.3 Potts model. 

Another important model which may described in terms of random curves in the Q-state 
Potts model. This is most easily considered on square lattice, at each site of which is a 
variable s{r) which can take Q (initially a positive integer) different values. The partition 
function is 

^PottS = Tr e^'^rr' Ur).ir') ^ Tv ^ (l - p + , (6) 

rr' 

with e^"^ — (1 — p)~^. The product may be expanded in a similar way to the case of the 
Ising model. All possible graphs Q will appear. Within each connected component of G 
the Potts spins must be equal, giving rise to a factor Q when the trace is performed. The 
result is _ 

^Potts = Ep'"(l-P)'''^3""'> (7) 

G 

where |^| is the number of edges in ^, \Q\ the number in its complement, and ||^|| is the 
number of connected components of G, which are called Fortuin-Kasteleyn (FK) clusters. 
This is the random cluster representation of the Potts model. When p is small, the mean 
cluster size is small. As p ^ Pc, it diverges, and for p > Pc there is an infinite cluster 
which contains a finite fraction of all the sites in the lattice. It should be noted that these 
FK clusters are not the same as the spin clusters within which the original Potts spins all 
take the same value. 

The limit Q ^ 1 gives another realisation of percolation - this time bond percolation 
on the square lattice. For Q ^ there is only one cluster. If at the same time x — > 
suitably, all loops are suppressed and the only graphs Q which contribute are spanning 
trees, which contain every site of the lattice. In the Potts partition function each possible 
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Figure 3: Example of FK clusters (heavy lines) in the random cluster representation of the 
Potts model, and the corresponding set of dense loops (medium heavy) on the medial lattice. 
The loops never cross the edges connecting sites in the same cluster. 

spanning tree is counted with the same weight, corresponding to the problem of uniform 
spanning trees (UST). The ensemble of paths on USTs connecting two points ri and r2 
turns out to be be that of loop-erased random walks. 

The random cluster model may be realised as a gas of dense loops in the way illustrated 
in Fig. ini These loops lie on the medial lattice, which is also square but has twice the 
number of sites. It may be shown that, at pc, the weights for the clusters are equivalent to 
counting each loop with a fugacity y/Q. Thus the boundaries of the critical FK clusters 
in the Q-state Potts model are the same in the scaling limit (if it exists) as the closed 
loops of the dense phase of the 0(n) model, with n = ^JQ. 

To generate an open path in the random cluster model connecting sites ri and r2 on the 
boundary we must choose 'wired' boundary conditions, in which p = 1 on all the edges 
parallel to the boundary, from r\ to r2, and free boundary conditions, with p = 0, along 
the remainder. 

2.4 Coulomb gas methods 

Many important results concerning the Q(n) model can be derived in a non-rigorous 
fashion using so-called Coulomb gas methods. For the purposes of comparison with later 
results from SLE, we now summarise these methods and collect a few relevant formulae. 
A much more complete discussion may be found in the review by Nienhuis|3]. 

We assume that the boundary conditions on the 0(n) spins are free, so that the partition 
function is a sum over closed loops only. First orient each loop at random. Rather than 
giving clockwise and anti-clockwise orientations the same weight n/2, give them complex 
weights e^^*^, where n = e^^^ + e~^*^ = 2cos6x- These may be taken into account, on 
the honeycomb lattice, by assigning a weight e^*^ at each vertex where an oriented loop 
turns R (resp. L). This transforms the non-local factors of n into local (albeit complex) 
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weights depending only on the local configuration at each vertex. 

Next transform to the height variables described above. By convention, the heights are 
taken to be integer multiples of tt. The local weights at each vertex now depend only 
on the differences of the three adjacent heights. The crucial assumption of the Coulomb 
gas approach is that, under the RG, this model fiows to one in which the lattice can be 
replaced by a continuum, and the heights go over into a gaussian free field, with partition 
function Z = J e"'^^'^' [dh] , where 

S=ig/4n) jiVkf^T. (8) 

As it stands, this is a simple free field theory. The height fiuctuations grow logarithmically: 
((/i(ri) — h{r2))^) ~ (2/(7) ln|ri — r2|, and the correlators of exponentials of the height 
decay with power laws: 

where Xq = /2g. All the subtleties come from the combined effects of the phase factors 
and the boundaries or the topology. This is particularly easy to see if we consider the 
model on a cylinder of circumference ^ and length L ^ i. In the simple gaussian model 
(jH} the correlation function between two points a distance L apart along the cylinder 
decays as exp ( — 2'nXqL/t). However, if x 7^ 0, loops which wrap around the cylinder are 
not counted correctly by the above prescription, because the total number of left turns 
minus right turns is then zero. We may arrange the correct factors by inserting e^''^*^'*/'^ 
at either end of the cylinder. This has the effect of modifying the partition function: one 
finds InZ ~ (7rc/6)(L/£) with 

c=l-6M!L)!. (10) 
9 

This dependence of the partition function is one way of determining the so-called central 
charge of the corresponding CFT (Sec. E)). The charges at each end of the cylinder also 
modify the scaling dimension Xq to {l/2g){{q — Qix/t^Y ~ (6^x/^)^)- 

The value of g may be fixed 17 in terms of the original discreteness of the height variables 
as follows: adding a term —A / cos2h d'^r to 5 in © ensures that, in the limit A —>■ 00, h 
will be an integer multiple of vr. For this deformation not to affect the critical behaviour, 
it must be marginal in the RG sense, which means that it must have scaling dimension 
X2 = 2. This condition then determines g = 1 — Ox/vr. 

2.4.1 Winding angle distribution 

A simple property which can be inferred from the Coulomb gas formulation is the winding 
angle distribution. Consider a cylinder of circumference 2n and a path that winds around 
it. What the probability that it winds through an angle 6 around the cylinder while 
it moves a distance L ^ 1 along the axis? This will correspond to a height difference 
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Ah = 71(6/271) between the ends of the cyhnder, and therefore an additional free energy 
{g/4:7T)(27iL){6/2Ly. The probabihty density is therefore 

P{e) oc expi^gOySL) , (11) 

so that 6 is normally distributed with variance {4/g)L. This result will be useful later 
(Sec. I3.6|l for comparison with SLE. 

2.4.2 A^-leg exponent 

As a final simple exponent prediction, consider the correlation function ($Ar(ri)$7v(r2)) 
of the iV-leg operator, which in the language of the 0(n) model is $Ar = s^^ ■ ■ ■ s^^, where 
none of the indices are equal. It gives the probability that mutually non-intersecting 
curves connect the two points. Taking them a distance L ^ i apart along the cylinder, 
we can choose to orient them all in the same sense, corresponding to a discontinuity in h 
of Ntt in going around the cylinder. Thus we can write h = ttNv/I + h, where < w < £ 
is the coordinate around the cylinder, and h{v + £) = h{v) . This gives 

(<l>;v(ri)$^(r2)) ~ e-(3/^")(^"/^)'^+('^^/'^^)-("^^/^^) . (12) 

The second term in the exponent comes from the integral over the fluctuations h, and 
the last from the partition function. They differ because in the numerator, once there 
are curves spanning the length of the cylinder, loops around it, which give the correction 
term in (fTUI) . are forbidden. (fT^ then gives 



XN = igNy8)~{g-lf/2g. 



(13) 
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3 



SLE 



3.1 



The postulates of SLE. 



SLE gives a description of the continuum limit of the lattice curves connecting two points 
on the boundary of a domain T> which were introduced in Sec. |2l The idea is to define a 
measure /x(7; V, ri, r2)) on these continuous curves. (Note that the notion of a probability 
density of such objects does not make sense, but the more general concept of a measure 
does.) 

There are two basic properties of this continuum limit which must either be assumed, or, 
better, proven to hold for a particular lattice model. The first is the continuum version 
of Property 1: 

Property 3.1 (continuum version). Denote the curve by 7, and divide it into two disjoint 
parts: 71 from ri to r, and 72 from r to r2. Then the conditional measure /i(72|7i; 1^, ri, r2) 
is the same as /i(72; T> \ 71, r, r2). 

This property we expect to be true for the scaling limit of all such curves in the 0{n) 
model (at least for n > 0), even away from the critical point. However the second property 
encodes the notion of conformal invariance, and it should be valid, if at all, only at x = Xc 
and, separately, for x > Xc- 

Property 3.2 (conformal invariance.) Let $ be a conformal mapping of the interior of 
the domain V onto the interior of T>', so that the points (ri, r2) on the boundary of T> are 
mapped to points (^"1,^2) on the boundary of V. The measure fi on curves in T> induces 
a measure $ * /i on the image curves in T>'. The conformal invariance property states that 
this is the same as the measure which would be obtained as the continuum limit of lattice 
curves from r[ to r'2 in V. That is 



3.2 Loewner's equation 

We have seen that, on the lattice, the curves 7 may be 'grown' through a discrete ex- 
ploration process. The Loewner process is the continuum version of this. Because of 
Property 2 it suffices to describe this in a standard domain V, which is taken to be the 
upper half plane H, with the points ri and r2 being the origin and infinity respectively. 

The first thing to notice is that, although on the honeycomb lattice the growing path 
does not intersect itself, in the continuum limit it might (although it still should not cross 
itself.) This means that there may be regions enclosed by the path which are not on 
the path but nevertheless are not reachable from infinity without crossing it. We call 
the union of the set of such points, together with the curve itself, up to time t, the hull 
Kt- (This is a slightly different usage of this term from that in the physics percolation 



($ * /i)(7; V, n, rs) = /i($(7); V\ r^, r^) . 



(14) 
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Figure 4: Schematic view of a trace and its hull. 



literature.) It is the complement of the connected component of the half plane which 
includes oo, itself denoted by H \ Kf. See Fig. HI 

Another property which often holds in the half-plane is that of reflection invariance: the 
distribution of lattice paths starting from the origin and ending at oo is invariant under 
X —X. For the lattice paths in the 0{n) model discussed in Sec. 12.21 this follows from 
the symmetry of the underlying weights, but for the boundaries of the FK clusters in the 
Potts model it is a consequence of duality. Not all simple curves in lattice models have 
this property. For example if we consider the 3-state Potts model in which the spins on 
the negative and positive real axes are fixed to different values, there is a simple lattice 
curve which forms the outer boundary of the spin cluster containing the positive real axis. 
This is not the same as the boundary of the spin cluster containing the negative real axis, 
and it is not in general symmetric under reflections. 

Since H \ i^^j is simply connected, by the Riemann mapping theorem it can be mapped 
into the standard domain H by an analytic function gt{z). Because this preserves the 
real axis outside Kt it is in fact real analytic. It is not unique, but can be made so by 
imposing the behaviour as z —>■ oo 



It can be shown that, as the path grows, the coefficient of 1/z is monotonic increasing 
(essentially it is the electric dipole moment of Kt and its mirror image in the real axis.) 
Therefore we may reparametrise time so that this coefficient is 2t. (The factor 2 is 
conventional.) Note that the length of the curve is not be a useful parametrisation in the 
continuum limit, since the curve is a fractal. 

The function gt{z) maps the whole boundary of Kf onto part of the real axis. In particular, 
it maps the growing tip Tt to a real point at. Any point on the real axis outside Kt remains 
on the real axis. As the path grows, the point at moves on the real axis. For the path to 
describe a curve, it must always grow only at its tip, and this means that the function at 
must be continuous, but not necessarily differentiable. 

A simple but instructive example is when 7 is a straight line growing vertically upwards 




(15) 
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from a fixed point a. In this case 

gt{z) = a + {{z - af + 4tf\ (16) 

This satisfies fjlSj) . and Tt = 2%\ft. More comphcated deterministic examples can be 
found jTH|. In particular, at oc t^/^ describes a straight line growing at a fixed angle to the 
real axis. 

Loewner's idea was to describe the path 7 and the evolution of the tip Tt in terms of 
the evolution of the conformal mapping gti^z). It turns out that the equation of motion 
for gt{z) is simple: 

dt gt{z) - at ' 

This is Loewner's equation. The idea of the proof is straightforward. Imagine evolving 
the path for a time t, and then for a further short time 6t. The image of Kt+st under gt 
is a short vertical line above the point at on the real axis. Thus we can write, using (jl6|) 

9t+sti^) -at + iigtiz) - a,)' + . (18) 

Differentiating with respect to 6t and then letting 6t —>■ 0, we obtain (fTTj) . 

Note that, even if at is not differentiable (as is the case for SLE), (fT7j) gives for each 
point zq a solution gt{zo) which is differentiable with respect to t, up to the time when 
gt{zo) = cit- This is the time when zq is first included in Kt- However, it is sometimes (see 
Sec. El) useful to normalise the Loewner mapping differently, defining gt{z) = gt{z) — at, 
which always maps the growing tip Tt to the origin. If at is not differentiable, neither is gt, 
and the Loewner equation should be written in differential form as dgt = {2dt/gt) — dat- 

Given a growing path, we can determine the hull Kt and hence, in principle, the function 
gt{z) and thereby at = gtijt)- Conversely, given at we can integrate (fTTj) to find gt{z) and 
hence in determine the curve (although proving that this inverse problem gives a curve is 
not easy.) 



3.3 Schramm-Loewner Evolution 

In the case that we are interested in, 7 is a random curve, so that at is a random continuous 
function. What is the measure on Of? This is answered by the following remarkable result, 
due to Schramm!^: 

Theorem. If Properties 3.1 and 3.2 hold, together with reflection symmetry, then at is 
proportional to a standard Brownian motion. 

That is 

at = y^Bt, (19) 

so that (at) = 0, ((aj^ — atj^) = K\ti — t2\. The only undetermined parameter is k, 
the diffusion constant. It will turn out that different values of k, correspond to different 
universality classes of critical behaviour. 
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Figure 5: A hull evolved from ao for time ti, then to infinity. The measure on the image of the 
rest of the curve under gt^ is the same, according to the postulates of SLE, as a hull evolved 
from at^ to oo. 

The idea behind the proof is once again simple. As before, consider growing the curve 
for a time ti, giving 71, and denote the remainder 7 \ 71 = 72- Property 3.1 tells us that 
the conditional measure on 72 given 71 is the same as the measure on 72 in the domain 
H \ Ktj^, which, by Property 3.2, induces the same measure on gti{l2) in the domain H, 
shifted by (see Fig. E}. In terms of the function at this means that the probability 
law of at — at^i for t > ti, is the same as the law of at-t^- This implies that all the 
increments A„ = a(^n+i)5t — o-nSt are independent identically distributed random variables, 
for all 5t > 0. The only process that satisfies this is Brownian motion with a possible 
drift term: at = ^/K,Bt + at. Reflection symmetry then implies that a = 0. 

3.4 Simple Properties of SLE 
3.4.1 Phases of SLE 

Many of the results discussed in this section have been proved by Rohde and Schramm|19j. 
First we address the question of how the trace (the trajectory of Tt) looks for different 
values of k. For k = 0, it is a vertical straight line. As k increases, the trace should 
randomly turn to the L or R more frequently. However, it turns out that there are 
quahtative differences at critical values of k. To see this, let us first study the process on 
the real axis. Let Xt = gt{xo) — at be the distance between the image at time t of a point 
which starts at xq and the image at of the tip. It obeys the stochastic equation 



Physicists often write such an equation as x = {2/x) — rjt where rjt is 'white noise' of 
strength k. Of course this does not make sense since Xt is not differentiable. Such 
equations are always to be interpreted in the Tto sense', that is, as the limit as (5t — >• of 
the forward difference equation Xt+st ~ a;t + {25t/xt) + Jt"^^^ rjfdt' . 

f|20|l is known as the Bessel process. (If we set Rt = {D — l)^/^a;t/2 and = A/{D — 1) 
it describes the distance Rt from the origin of a Brownian particle in D dimensions.) 




(20) 
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Figure 6: The trace is about to hit the axis at xq and enclose a region. At the time this happens, 
the whole region including the point xq is mapped by gt to the same point at. 



The point Xt is repelled from the origin but it is also subject to a random force. Its 
ultimate fate can be inferred from the following crude argument: if we ignore the random 
force, xf ~ 4t, while, in the absence of the repulsive term, (xf) ~ nt. Thus for k < 4 
the repulsive force wins and the particle escapes to infinity, while for k > 4 the noise 
dominates and the particle collides with the origin in finite time (at which point the 
equation breaks down.) A more careful analysis confirms this. What does this collision 
signify in terms of the behaviour of the trace? In Fig. IHl we show a trace which is about 
to hit the real axis at the point xq, thus engulfing a whole region. This is visible from 
infinity only through a very small opening, which means that, under gt, it gets mapped 
to a very small region. In fact, as the tip Tt approaches xq, the size of the image of this 
region shrinks to zero. When the gap closes, the whole region enclosed by the trace, as 
well as Tt and Xq, are mapped in to the single point at, which means, in particular, that 
Xt ^ 0. The above argument shows that for /t < 4 this never happens: the trace never 
hits the real axis (with probability 1.) For the same reason, it neither hits itself. Thus 
for K < 4 the trace 7 is a simple curve. 

The opposite is true for k > 4: points on the real axis are continually colliding with the 
image at of the tip. This means that the trace is continually touching both itself and the 
real axis, at the same time engulfing whole regions. Moreover, since it is self-similar object, 
it does this on all scales, an infinite number of times within any finite neighbourhood! 
Eventually the trace swallows the whole half plane: every point is ultimately mapped into 
at. For K < 4 only the points on the trace itself suffer this fate. The case k = 4 is more 
tricky: in fact the trace is then also a simple curve. 

When K is just above 4, the images of points on the real axis under gt collide with at only 
when there happen to be rare events when the random force is strong enough to overcome 
the repulsion. When this happens, whole segments of the real axis are swallowed at one 
time, corresponding to the picture described above. Conversely, for large k, the repulsive 
force is negligible except for very small Xt. In that case, two different starting points move 
with synchronised Brownian motions until the one which started off closer to the origin is 
swallowed. Thus the real line is eaten up in a continuous fashion rather than piecemeal. 
There are no finite regions swallowed by the trace which are not on the trace itself. This 
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means that the trace is space-filling: 7 intersects every neighbourhood of every point in 
the upper half plane. We shall argue later (Sec. 14.3. Ij) that the fractal dimension of the 
trace is df = 1 + k,/8 for k < 8 and 2 for k > 8. Thus it becomes space-filling for all 

K>8. 

3.4.2 SLE duality 

For K > 4 the curve continually touches itself and therefore its hull Kf contains earlier 
portions of the trace (see Fig.|3]). However, the frontier of Kf (i.e. the boundary of H\i^f, 
minus any portions of the real axis), is by definition a simple curve. A beautiful result, 
first suggested by Duplantierj^U], and proved by Beffara'^ for the case k = 6, is that 
locally this curve is an SLE^j, with 

k = 16/K. (21) 

For example, the exterior of a percolation cluster contains many 'fjords' which, on the 
lattice, are connected to the main ocean by a neck of water which is only a finite number of 
lattice spacings wide. These are sufficiently frequent and the fjords macroscopically large 
that they survive in the continuum limit. SLEg describes the boundaries of the clusters, 
including the coastline of all the fjords. However, the coastline as seen from the ocean is 
a simple curve, which is locally SLEg/s, the same as that conjectured for a self-avoiding 
walk. This suggests, for example, that locally the frontier of a percolation cluster and a 
self-avoiding walk are the same in the scaling limit. In Sec. El we show that SLE^ and 
SLEk correspond to CFTs with the same value of the central charge c. 

3.5 Special values of k 
3.5.1 Locality 

[This subsection and the next are more technical and may be omitted at a first reading.] 
We have defined SLE in terms of curves which connect the origin and infinity in the upper 
half plane. Property 2 then allows us to define it for any pair of boundary points in any 
simply connected domain, by a conformal mapping. It is interesting to study how the 
variation of the domain affects the SLE equation. Let A be a simply connected region 
connected to the real axis which is at some finite distance from the origin (see Fig. [7j). 
Consider a trace 74, with hull Kt, which grows from the origin according to SLE in the 
domain H \ A. According to Property 2, we can do this by first making a conformal 
mapping Hq which removes A, and then a map gt which removes the image Kf = hQ{Kt). 
This would be described by SLE in ho(H. \ A), except that the Loewner 'time' would not 
in general be the same as t. However, another way to think about this is to first use a 
SLE map gt in H to remove Kt, then another map, call it ht, to remove gt{A). Since both 
these procedures end up removing KtUA, and all the maps are assumed to be normalised 
at infinity in the standard way (jl5p . they must be identical, that is ht o g^ = g^ o h^ (see 
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at at 

Figure 7: An SLE hull in H \ A and two different ways of removing it: either by first removing 
A through Hq and then using a Loewner map gt in the image of H \ A; or by removing Kt first 
with gt and then removing the image of A with hf. Since all maps are normalised, this diagram 
commutes. 



Fig.[7j). If gt maps the growing tip Tt to at, then after both mappings it goes to dt = ht{at). 
We would like to understand the law of dt- 

Rather than working this out in full generality (see for example jTHI), let us suppose that 
A is a short vertical segment {x,x + ie) with e <^ x, and that t = dtis infinitesimal. Then, 
under g^t, x ^ x + 2dt/x and e — e(l — 2dt/x'^). The map that removes this is (see p6|) ') 

hat{z) = ((z - X - 2dt/xY + £^(1 - 2dt/x^f) ^'"^ + x + 2dt/x . (22) 

In order to find ddt-, we need to set z = adt = \fi^dBt in this expression. Carefully 
expanding this to first order in dt, remembering that {dBtY = dt, and also taking the 
first non-zero contribution in e/x, gives after a few lines of algebra 

ddt = (1 - e^/x^)^dBt + Q){e^/x^)dt . (23) 

The factor in front of the stochastic term may be removed by rescaling dt: this restores 
the correct Loewner time. But there is also a drift term, corresponding to the effect of A. 
For K < 6 we see that the SLE is initially repelled from A. From the point of view of the 
exploration process for the Ising model discussed in Sec. I2.1.H this makes sense: if the 
spins along the positive real axis and on A are fixed to be up, then the spin just above 
the origin is more likely to be up than down, and so 7 is more likely to turn to the left. 

For K = 6, however, this is no longer the case: the presence of A does not affect the initial 
behaviour of the curve. This is a particular case of the property of locality when k = 6, 
which states that, for any A as defined above, the law of in H \ A is, up to a time 
reparametrisation, the same as the law of Kt in H, as long as KtCiA = 0. That is, up to the 
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time that the curve hits A, it doesn't know it's there. Such a property would be expected 
for the cluster boundaries of uncorrelated Ising spins on the lattice, i.e. percolation. This 
is then consistent with the identification of percolation cluster boundaries with SLEg. 

3.5.2 Restriction 

It is also interesting to work out how the local scale transforms in going from at to fij. A 
measure of this is h[{at). A similar calculation starting from ()22|1 gives, in the same limit 
as above, 

d{h\at)) = h'M - h'oiO) = {e^/x')^dBt + \ ({e'/x') + {k - |)(3eVa;^)) . (24) 

Now something special happens when /t = |. The drift term in d{h'{at)) does not then 
vanish, but if we take the appropriate power d(h[{atY^^) does. This implies that the 
mean of h[{atY^^ is conserved. Now at t = it takes the value $^(0)^/*, where = 
ho is the map that removes A. If Kt hits A at time T it can be seen from (|22j) that 
liuit^T ht{atY^^ = 0. On the other hand, if it never hits A then limj^oo ^K^t)^'''^ = 1- 
Therefore $^(0)^/^ gives the probability that the curve 7 does not intersect A. 

This is a remarkable result in that it depends only on the value of $^ at the starting point 
of the SLE (assuming of course that $a is correctly normalised at infinity.) However it has 
the following even more interesting consequence. Let ^a{z) = ^a{z) — ^{0). Consider the 
ensemble of all SLEg/s in H, and the sub-ensemble consisting of all those curves 7 which 
do not hit A. Then the measure on the image $a(7) in H is again given by SLEg/a. The 
way to show this is to realise that the measure on 7 is characterised by the probability 
P(7 n A' = 0) that 7 does not hit A' for all possible A'. The probability that $^(7) does 
not hit A', given that 7 does not hit A, is the ratio of the probabilities P{'jr]^A^{A') = 0) 
and P(7 fl A = 0). By the above result, the first factor is the derivative at the origin of 
the map $a' ° which removes A then A', while the second is the derivative of the map 
which removes A. Thus 

5/8 

p{^a{i) n A' = 0|7 n a = 0) = f (^a^^IaIM] = $^,(0)5/8 = p(7 n A' = 0) . (25) 

V $a(o) / 

Since this is true for all A', it follows that the law of ^"^(7) given that 7 does not intersect A 
is the same as that of 7. This is called the restriction property. Note that while, according 
to Property 2, the law of an SLE in any simply connected subset of H is determined by 
the conformal mapping of this subset to H, the restriction property is stronger than this, 
and it holds only when = ^■ 

We expect such a property to hold for the continuum limit of self-avoiding walks, assuming 
it exists. On the lattice, every walk of the same length is counted with the same weight. 
That is, the measure is uniform. If we consider the sub-ensemble of such walks which 
avoid a region A, the measure on the remainder should still be uniform. This will be true 
if the restriction property holds. This supports the identification of self-avoiding walks 
with SLEg/3. 



SLE for Theoretical Physicists 



22 



3.6 Radial SLE and the winding angle 

So far we have discussed a version of SLE that gives a conformally invariant measure on 
curves which connect two distinct boundary points of a simply connected domain V. For 
this reason it is called chordal SLE. There are variants which describe other situations. 
For example, one could consider curves 7 which connect a given point ri on the boundary 
to an interior point r2. The Riemann mapping theorem allows us to map conformally onto 
another simple connected domain, with r2 being mapped to any preassigned interior point. 
It is simplest to choose for the standard domain the unit disc U, with r2 being mapped to 
the origin. So we are considering curves 7 which connect a given point e*^° on the boundary 
with the origin. As before, we may consider growing the curve dynamically. Let Kt be 
the hull of that portion which exists up to time t. Then there exists a conformal mapping 
Qt which takes \J \Kt to U, such that gt{0) = 0. There is one more free parameter, which 
corresponds to a global rotation: we use this to impose the condition that ^'^(O) is real 
and positive. One can then show that, as the curve grows, this quantity is monotonically 
increasing, and we can use this to reparametrise time so that (7^(0) = e*. This normalised 
mapping then takes the growing tip Tt to a point e*^' on the boundary. 

Loewner's theorem tells us that gt{z) / gt{z) , when expressed as a function of gt{z), should 
be holomorphic in U apart from a simple pole at e*^'. Since gt preserves the unit circle 
outside Kt, gt{z)/gt{z) should be pure imaginary when = 1, and in order that 

g'tiO) = e*, it should approach 1 as gt{z) 0. The only possibility is 

gt{z)+e''^ 

~ir - gt{z) - ■ ^^^^ 

This is the radial Loewner equation. In fact this is the version considered by LoewnerjJ]. 

It can now be argued, as before, that given Properties 1 & 2 (suitably reworded to cover 
the case when r2 is an interior point) together with reflection, Qt must be proportional to 
a standard Brownian motion. This defines radial SLE. It is not immediately obvious how 
the radial and chordal versions are related. However, it can be shown that, if the trace 
of radial SLE hits the boundary of the unit disc at e*^*i at time ti, then the law of Kt in 
radial SLE, for t < ti, is the same chordal SLE conditioned to begin at e*^*^"-* and end at 
e*^'i, up to a reparametrisation of time. This assures us that, in using the chordal and 
radial versions with the same we are describing the same physical problem. 

However, one feature that the trace of radial SLE possesses which chordal SLE does not 
is the property that it can wind around the origin. The winding angle at time t is simply 
9t — 9Q. Therefore it is normally distributed with variance nt. At this point we can make 
a connection to the Coulomb gas analysis of the 0{n) model in Sec. 12.4.11 where it was 
shown that the variance in the winding angle on a cylinder of length L is asymptotically 
{4/g)L. A semi-infinite cylinder, parametrised by w, is conformally equivalent to the unit 
disc by the mapping z = e"*". Asymptotically, Hew lie w — t under Loewner evolution. 
Thus we can identify L t and hence 

K = A/g. (27) 
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3.6.1 Identification with lattice models 

This result allows use to make a tentative identification with the various realisations of the 
0(n) model described in Sec. 12.21 We have, using (j77j) . n = —2 cos(47r/fi;) with 2 < k < 4 
describing the critical point at Xc, and 4 < n < 8 corresponding to the dense phase. Some 
important special cases are therefore: 

• K = —2: loop-erased random walks (proven in |24j): 

• K = |: self-avoiding walks, as already suggested by the restriction property. Sec. 13. 5. 2] 
unproven, but see for many consequences; 

• K = 3: cluster boundaries in the Ising model, however as yet unproven; 

• K = 4: BCSOS model of roughening transition (equivalent to the 4-state Potts 
model and the double dimer model), as yet unproven; also certain level lines of a 
gaussian random field and the 'harmonic explorer' (proven in [211); also believed to 
be dual to the Kosterlitz-Thouless transition in the XY model; 

• n = 6: cluster boundaries in percolation (proven in 0); 

• K = 8: dense phase of self-avoiding walks; boundaries of uniform spanning trees 
(proven in ^^). 

It should be noted that no lattice candidates for k > 8, or for the dual values k, < 2, have 
been proposed. This possibly has to do with the fact that, for k, > 8, the SLE trace is 
not reversible: the law on curves from ri to r2 is not the same as the law obtained by 
interchanging the points. Evidently, curves in equilibrium lattice models should satisfy 
reversibility. 
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4 Calculating with SLE 



SLE shows that the measure on the continuum hmit of single curves in various lattice 
models is given in terms of Id Brownian motion. However, it is not at all clear how 
thereby to deduce interesting physical consequences. We first describe two relatively 
simple computations in two-dimensional percolation which can be done using SLE. 

4.1 Schramm's formula 

Given a curve 7 connecting two points ri and r2 on the boundary of a domain V, what is 
the probability that it passes to the left of a given interior point? This is not a question 
which is natural in conventional approaches to critical behaviour, but which is very simply 
answered within SLE|25j. 

As usual, we can consider V to be the upper half plane, and take ri = and r2 to be 
at infinity. The curve is then described by chordal SLE starting at a^. Label the interior 
point by the complex number (. 

Denote the probability that 7 passes to the left of ( by P{C, C; cto) (we include the depen- 
dence on ( to emphasise the fact that this is a not a holomorphic function.) Consider 
evolving the SLE for an infinitesimal time dt. The function gdt will map the remainder of 
7 into its image 7', which, however, by Properties 1 & 2, will have the same measure as 
SLE started from a^t = + y/RdBt. At the same time, C 9dt{C) = C + '^dt/{( — oq). 
Moreover, 7' lies to the left of (' iff 7 lies to the left of (. Therefore 



where the average (. . .) is over all realisations of Brownian motion dBf up to time dt. 
Taylor expanding, using {dBt) = and {{dBtY) = dt, and equating the coefficient of dt 
to zero gives 



probabilities in stochastic differential equations. 

By scale invariance P in fact depends only on the angle 9 subtended between ( ~ (^0 and 
the real axis. Thus ()29|) reduces to an ordinary second-order linear differential equation, 
which is in fact hypergeometric. The boundary conditions are that P = and 1 when 
6 = IT and respectively, which gives (specialising to k = 6) 



P(C, C; «o) = {P{C + 2dt/{C - ao), C + 2rft/(C - ao), ao + V^dB^)) , 



(28) 




P 




(cot^)2Fi( 



12 3. 

2' 3' 2' 



COt^ 9) . 



(30) 



Note that this may also be written in terms of a single quadrature since one solution of 
is P = const. 
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Figure 8: Is there a crossing on the white discs from (0,X2) to (— oo,xi)? This happens if and 
only if xi gets swallowed by the SLE before X2 ■ 

4.2 Crossing probability 

Given a critical percolation problem inside a simply connected domain V, what is the 
probability that a cluster connects two disjoint segments AB and CD of the boundary? 
This problem was conjectured to be conformally invariant and (probably) first studied 
numerically in [2ni- A formula based on CFT as well as a certain amount of guesswork 
was conjectured in It was proved, for the continuum limit of site percolation on the 
triangular lattice, by Smirnov[7j. 

Within SLE, it takes a certain amount of ingenuity(S] to relate this problem to a question 
about a single curve. As usual, let V be the upper half plane. It is always possible 
to make a fractional linear conformal mapping which takes AB into (— oo,a;i) and CD 
into (0,X2), where xi < and X2 > 0. Now go back to the lattice picture and consider 
critical site percolation on the triangular lattice in the upper half plane, so that each 
site is independently coloured black or white with equal probabilities |. Choose all the 
boundary sites on the positive real axis to be white, all those on the negative real axis 
to be black (see Fig. |H1). There is a cluster boundary starting at the origin, which, in the 
continuum limit, will be described by SLEg. Since k > 4, it repeatedly hits the real axis, 
both to the L and R of the origin. Eventually every point on the real axis is swallowed. 
Either xi is swallowed before X2, or vice versa. 

Note that every site on the L of the curve is black, and every site on its R is white. 
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Suppose that xi is swallowed before X2- Then, at the moment it is swallowed, there exists 
a continuous path on the white sites, just to the R of the curve, which connects (0, X2) to 
the row just above (—00, xi). On the other hand, if X2 is swallowed before xi, there exists 
a continuous path on the black sites, just to the L of the curve, connecting 0— to a point 
on the real axis to the R of 0:2. This path forms a barrier (as in the game of Hex) to the 
possibility of a white crossing from (0,0:2) to (— oo,a;i). Hence there is such a crossing if 
and only if xi is swallowed before X2 by the SLE curve. 

Recall that in Sec. 13.4.11 we related the swallowing of a point Xq on the real axis to the 
vanishing of Xt = gt{xt) — o-t, which undergoes a Bessel process on the real line. Therefore 

Pr(crossing from (0,^2) to (— oo,Xi)) = Pr(xit vanishes before X2t) ■ (31) 

Denote this by P{xi,X2). By generalising the SLE to start at rather than 0, we can 
write a differential equation for this in similar manner to (j29p : 

(^—^ + + Sl^l ^(^1' ao) = . (32) 

\Xi—aQOXi X2 — aoOX2 2 oaQ J 

Translational invariance implies that we can replace da^ by —{dx^ + 0^2) ■ Finally, P can 
in fact depend only on the ratio rj = {x2 — ao) /(oq — Xi), which again reduces the equation 
to hypergeometric form. The solution is (specialising to k = 6 for percolation) 

^ = f(P|Ty^^^i^i(ii4;^)- (33) 

It should be mentioned that this is but one of a number of percolation crossing formulae. 
Another, conjectured by Watts }27j. for the probability that there is cluster which simul- 
taneously connects AB to CD and BC to DA, has since been proved by Dubedat j28j . 
However, other formulae, for example for the mean number of distinct clusters connect- 
ing AB and CD [211, and for the probability that exactly N distinct clusters cross an 
annuluspU]. yet unproven using SLE methods. 



4.3 Critical exponents from SLE 

Many of the critical exponents which have previously been conjectured by Coulomb gas 
or CFT methods may be derived rigorously using SLE, once the underlying postulates 
are assumed or proved. However SLE describes the measure on just a single curve, and 
in the papers of LSW a great deal of ingenuity has gone into showing how to relate this 
to all the other exponents. There is not space in this article to do these justice. Instead 
we describe three examples which give the flavour of the arguments, which initially may 
appear quite unconventional compared with the more traditional approaches. 
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4.3.1 The fractal dimension of SLE 

The fractal dimension of any geometrical object embedded in the plane can be defined 
roughly as follows: let N{e) be the minimum number of small discs of radius e required 
to cover the object. Then if N{e) ~ e~'^f as e — > 0, df is the fractal dimension. 

One way of computing dj for a random curve 7 in the plane is to ask for the probability 
P{x, y, e) that a given point ( = x + iy lies within a distance e of 7. A simple scaling 
argument shows that if P behaves like e^f{x,y) as e — > 0, then 6 = 2 — df. We can 
derive a differential equation for P along the lines of the previous calculation. The only 
difference is that under the conformal mapping gat, e — >■ \g'dt{0\^ ~ (1 — 2dtRe{l/('^))e. 
The differential equation (written for convenience in cartesian coordinates) is 



Now P is dimensionless and therefore should have the form (e/r)^"'^^ times a function of 
the polar angle 6. In fact, the simple ansatz P = e^~'^^?/"(x^ + y^)^, with a + 2f3 = df — 2 
satisfies the equation. [The reason this works is connected with the simple form for the 
correlator ($202,i02,i) discussed in Sec. 15.4. fl ] This gives a = (k — 8)^/8k, j3 = {k — 8)/2k 
and 



This is correct for k < 8: otherwise there is another solution with a = (3 = and df = 2. 
A more careful statement and proof of this result can be found in |31j . 

We see that the fractal dimension increases steadily from the value 1 when k = (corre- 
sponding to a straight line) to a maximum value of 2 when k, = 8. Beyond this value 7 
becomes space-filling: every point in the upper half plane lies on the curve. 

4.3.2 Crossing exponent 

Consider a critical percolation problem in the upper half plane. What is the asymptotic 
behaviour as r ^ 00 of the probability that the interval (0, 1) on the real axis is connected 
to the interval (r, 00)? We expect this to decay as some power of r. The value of this 
exponent may be found by taking the appropriate limit of the crossing formula ()33p . but 
instead we shall compute it directly. In order for there to be a crossing cluster, there must 
be two cluster boundaries which also cross between the two intervals, and which bound 
this cluster above and below. Denote the upper boundary by 7. Then we need to know 
the probability P{r) of there being another spanning curve lying between 7 and {l,r), 
averaged over all realisations of 7. Because of the locality property, the measure on 7 
is independent of the existence of the lower boundary, and is given by SLEg conditioned 
not to hit the real axis along (1,t). Note that because k > 4 it will eventually hit the 
real axis at some point to the right of r. For this reason we can do the computation for 
general k > 4, although it gives the actual crossing exponent only if /t = 6. 




(34) 



df = 1 + k/8 . 



(35) 
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Consider the behaviour of P{r) under the conformal mapping gdt{z) ~ z+{2dt/ z) — y/ndBt 
(which maps the growing tip into 0.) The crossing probabihty should be conformally 
invariant and depend only on the ratio of the lengths of the two intervals, hence, by an 
argument which by now should be familiar, 

P(r) = {P{g,,{r)/g,,{l))) . (36) 

Expanding this out, remembering as usual that {dBtY = dt, and setting to zero the 0{dt) 
term, we find for r ^ 1 

{k - 2)rP'{r) + \Kr'^P"{r) = , (37) 
with the solution P(r) oc r^'^'^^^^/'^ for k> A. Setting k = Q then gives the result |. 

4.3.3 The one-arm exponent 

Consider critical lattice percolation inside some finite region (for example a disc of radius 
R). What is the probability that a given site (e.g. the origin) is connected to a finite 
segment 5* of the boundary? This should decay like R~^, where A is sometimes called 
the one-arm exponent. If we try to formulate this in the continuum, we immediately 
run up against the problem that all clusters are fractal with dimension < 2, and so the 
probability of any given point being in any given cluster is zero. Instead, one may ask 
about the probability P{r) that the cluster connected to S gets within a distance r of the 
origin. This should behave like (r/R)^. We can now set i? = 1 and treat the problem 
using radial SLEg. 

Consider now a radial SLE^ which starts at e*^". If k > 4 it will continually hit the 
boundary. Let P{9 — 9o,t) be the probability that the segment {0q,9) of the boundary 
has not been swallowed by time t. Then, by considering the evolution as usual under g^t, 

p{e, Oo, t) = {P{e + de, Oq + dOo, t - dt)) , (38) 

where dO = cot {{6 — 9Q)/2)dt and dO^ = ^/KdEf. Setting 6q = and equating to zero the 
0{dt) term, we find the time-dependent differential equation 

dtP = cot{e/2)deP + \ndlP . (39) 

This has the form of a backwards Fokker-Plank equation. 

Now, since g[{^) = e*, it is reasonable that, after time t, the SLE gets within a distance 
0(e~*) of the origin. Therefore we can interpret P as roughly the probability that the 
cluster connected to (0, 6) gets within a distance r ~ e~* of the origin. A more careful 
argument [^2] confirms this. The boundary conditions are P(0, t) = as — 0, and (with 
more difficulty) dgP{9,t) = aX 9 = 27t. The solution may then be found by inspection 
to be 

Poce~^*(sin(^/4))'-^/", (40) 
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where A = (k^ — 16)/32k. For percolation this gives ^, in agreement with Coulomb gas 
arguments 0. 

The appearance of differential operators such as that in ()39|) will become clear from the 
CFT perspective in Sec. I5.4.1I If instead of choosing Neuman boundary conditions at 
6 = 2n we impose P = 0, the same equation gives the bulk 2-leg exponent X2, which is 
also related to the fractal dimension hj df = 2 — X2. 



SLE for Theoretical Physicists 



30 



5 Relation to conformal field theory 
5.1 Basics of CFT. 

The reader who aheady knows a httle about CFT will have recognised the differential 
equations in Sec. E]as being very similar in form to the BPZ equations[2SI satisfied by the 
correlation functions of a 02, i operator, corresponding to a highest weight representation 
of the Virasoro algebra with a level 2 null state. 

For those readers for whom the above paragraph makes no sense, and in any case to make 
the argument self-contained, we first introduce the concepts of (boundary) conformal 
field theory (BCFT.) We stress that these are heuristic in nature - they serve only as 
a guide to formulating the basic principles of CFT which can then be developed into a 
mathematically consistent theory. For a longer introduction to BCFT see jSl] and, for a 
complete account of CFT, 35j. 

We have at the back of our minds a euclidean field theory defined as a path integral 
over some set of fundamental fields {i'{r)}. The partition function is Z = J e'^^^^^'^ldip] 
where the action S{{ip}) = J-jy C{{ilj})d'^ is an integral over a local lagrangian density. 
These fields may be thought of as smeared-out continuum versions of the lattice degrees of 
freedom. As in any field theory, this continuum limit involves renormalisation. There are 
so-called local scaling operators (pf^ (r) which are particular functionals of the fundamental 
degrees of freedom, which have the property that we can define renormalised scaling 
operators 4>j{r) = a~^^ (j)^^\r) whose correlators are finite in the continuum limit a 0, 
that is 

linia~E,-.(0(o)(^^) . . ^ ^^^(^^) _ ..MrN)) (41) 

exists. The numbers xj are called the scaling dimensions, and are related to the various 
critical exponents. They are related to the conformal weights (hj,hj) by Xj = hj + hj; 
the difference hj — hj = sj is called the spin of (pj, and describes its behaviour under 
rotations. There are also boundary operators, localised on the boundary, which have only 
a single conformal weight equal to their scaling dimension. 

The theory is developed independently of any particular set of fundamental fields or 
lagrangian. An important role in this is played by the stress tensor T^'^{r), defined as the 
local response of the action to a change in the metric: 

6S = {1/A7r) f T'^'Sg^.d^r. (42) 
Jv 

Invariance under local rotations and scale transformations usually implies that T'^'^ is 
symmetric and traceless: = 0. This also implies invariance under conformal transfor- 
mations, corresponding to 6g^u oc f{r)gfj,u. 

In two dimensional flat space, infinitesimal coordinate transformations r"^ = 

r'^ + a'^(r) correspond to infinitesimal transformations of the metric with Sg'^'^ = —{d^a'^ + 
d'^a"). It is important to realise that under these transformations the underlying lattice. 
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or UV cut-off, is not transformed. Otherwise they would amount to a trivial reparametri- 
sation. For a conformal transformation, a'^(r) is given by an analytic function: in complex 
coordinates {z,z), dzO.^ = 0, so = a{z) is holomorphic. However, such a function can- 
not be small everywhere (unless it is constant), so it is necessary to consider coordinate 
transformations which are not everywhere conformal. 

Consider therefore two concentric semicircles Fi and F2 in the upper half plane, centred 
on the origin, and of radii Ri < R2. For |r| < Ri let a'^ be conformal, with = a{z), 
while for \r\ > R2 take a'^ = 0. In between, is not conformal, but is differentiable, so 
that 5S = (— l/27r) /R^<|r|<_R2 T'^^d^aydPr. This can be integrated by parts to give a term 
(l/27r) jRj^<:\r\<:R2dfj.T'^'^au(fr (which must vanish because is arbitrary in this region, 
implying that d^T^^ = 0) and two surface terms. That on F2 vanishes because = 
there. We are left with 

6S = (l/27r) / T^^^af^e^^dh , (43) 
where di^ is the line element along Fi. 

The fact that T'^" is conserved means, in complex coordinates, that dzTzz = d^T^z = 0, so 
that the correlations functions of T{z) = Tzz are holomorphic functions of z, while those 
of T = Tzz are antiholomorphic. may then be written 

6S=(l/2m)f T(z)a(z)dz + c.c.. (44) 
-'ri 

In any field theory with a boundary, it is necessary to impose some boundary condition. 
It can be argued that any translationally invariant boundary condition flows under the 
RG to conditions satisfying Txy = 0, which in complex coordinates means that T = T on 
the real axis. This means that the correlators of T are those of T analytically continued 
into the lower half plane. The second term in ()44p may then be dropped if the contour in 
the first term is around a complete circle. 

The conclusion of all this is that the effect of an infinitesimal conformal transformation 
on any correlator of observables inside Fi is the same as inserting a contour integral 
J T{z)a{z)dz /27ci into the correlator. 

Another important element of CFT is the operator product expansion (OPE) of the stress 
tensor with other local operators. Since T is holomorphic, this has the form 

T(^) ■ 0(0) = 5: V-HO) , (45) 

n 

where the 0*-"'' are (possibly new) local operators. By taking a{z) oc z (corresponding to 
a scale transformation) we see that (p^^^ = hep, where h is its scaling dimension. Similarly, 
by taking a = const., (f)^~^^ = dx4>- Local operators for which 0*-"'' vanishes for n > 1 are 
called primary. T itself is not primary: its OPE with itself takes the form 

T{z) ■ T(0) = c/2z' + {2/z')T{0) + (l/z)9,r(0) + ■ ■ ■ , (46) 

where c is the conformal anomaly number, ubiquitous in CFT. For example, the partition 
function on a long cylinder of length L and circumference i behaves as exp(7rcL/£), cf. (fTTl|) . 
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5.2 Radial quantisation 

This is the most important concept in understanding the hnk between SLE and CFT. 
We introduce it in the context of boundary CFT. As before, suppose there is some set of 
fundamental fields {'(p{r)}, with a Gibbs measure e~^^'^^[d-ip]. Let F be a semicircle in the 
upper half plane, centered on the origin. The Hilbert space of the BCFT is the function 
space (with a suitable norm) of field configurations {ipr} on F. 

The vacuum state is given by weighting each state lip'^) by the (normalised) path integral 
restricted to the interior of F and conditioned to take the specified values ip'^ on the 
boundary: 

|0) = / Wr] I , m e-^''^^ Wv) ■ (47) 

Note that because of scale invariance different choices of the radius of F are equivalent, 
up to a normalisation factor. 

Similarly, inserting a local operator (/)(0) at the origin into the path integral defines a 
state 10). This is called the operator-state correspondence of CFT. If we also insert 
(l/27ri) z^^^Tiz^dz^ where C lies inside F, we get a state The L„ act linearly on 

the Hilbert space. From the OPE we see that oc l^*-"-*), and that, in particular, 

Lo|0) = ^</>|0). If is primary, = for n > 1. From the OPE pHjl of T with itself 

follow the commutation relations for the L„ 

Lm] = (n - m)Ln+m + 3^072(72^ - l)Sn-m , (48) 

which are known as the Virasoro algebra. The state |0) together with all its descendants, 
formed by acting on |0) an arbitrary number of times with the L„ with n < — 1, give a 
highest weight representation (where the weight is defined as the eigenvalue of —Lq.) 

There is another way of generating such a highest weight representation. Suppose the 
boundary conditions on the negative and positive real axes are both conformal, that is 
they satisfy T = T, but they are different. The vacuum with these boundary conditions 
gives a highest weight state which it is sometimes useful to think of as corresponding to 
the insertion of a 'boundary condition changing' (bcc) operator at the origin. An example 
is the continuum limit of an Ising model in which the spins on the negative real axis are 
— 1, and those on the positive axis are +1. 

5.3 Curves and states. 

In this section we describe a way of associating states in the Hilbert space of the BCFT 
with the growing curves of the Loewner process. This was first understood by M. Bauer 
and D. Bernard[10j, but we shall present the arguments slightly differently. 

The boundary conditions associated with a bcc operator guarantee the existence, on the 
lattice, of a domain wall connecting the origin to infinity. Given a particular realisation 
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7, we can condition the Ising spins on its existence. We would like to be able to assume 
that this property continues to hold in the continuum limit: that is, we can condition the 
fields {ip} on the existence of a such a curve. However, this involves conditioning on an 
event with probability zero: it turns out that in general the probability that, with respect 
to the measure in the path integral, the probability that a domain wall hits the real axis 
somewhere in an interval of length e vanishes like e^. In what follows we shall regard e as 
small but fixed, and assume that the usual properties of SLE are applicable to this more 
general case. 

Any such curve may be generated by a Loewner process: denote as before the part of 
the curve up to time t by The existence of this curve depends on only the field 
configurations ijj in the interior of F, as long as 7t lies wholly inside this region. Then 
we can condition the fields contributing to the path integral on the existence of •jt, thus 
defining a state 

\lt)=IWT\l me-'^^^W^). (49) 

The path integral (over the whole of the upper half plane, not just the interior of F), when 
conditioned on 7^, gives a measure d^{^t) on these curves. The state 

\h) = \ht)^ Jdf,{^t)\lt) (50) 

is in fact independent of t, since it is just given by the path integral conditioned on there 
being a curve connecting the origin to infinity, as guaranteed by the boundary conditions. 
In fact we see that \h) is just the state corresponding to a boundary condition changing 
operator at the origin. 

However, d^{^t) is also given by the measure on at in Loewner evolution, through the 
iterated sequence of conformal mappings satisfying dgt = 2dt/gt — dat- This corresponds 
to an infinitesimal conformal mapping of the upper half plane minus Kt. As explained 
in the previous section, dgt corresponds to inserting {l/2m) J^{2dt/z — dat)T{z)dz. In 
operator language, this corresponds to acting on |7f) with 2L_2dt — L^idat where L„ = 
(l/27ri) Jcz''+^T{z)dz. Thus, for any < t, 

\gtAlt)) = Texp (^^*' i2L_2dt' - L.^da,)) \^t) , (51) 
where T denotes a time-ordered exponential. 

The measure on 74 is the product of the measure of 7f \ 7^1 , conditioned on 7f ^ , with the 
unconditioned measure on 7^^. The first is the same as the unconditioned measure on 
gtiilt), and the second is given by the measure on af for t' G [0,ti]. Thus we can rewrite 
both the measure and the state in as 

\ht) = J d^,{gtAlt)) I ^^/i(a,.t'e[OA])Te^ (2^-^*'"^-'^'^'') |^,^(^^)) . (52) 

For SLE, at is proportional to a Brownian process. The integration over realisations of this 
for t' G [0, ti] may be performed by breaking up the time interval into small segments of 



SLE for Theoretical Physicists 



34 



size St, expanding out the exponential to 0{St), using {BstY ~ St, and re-exponentiating. 
The result is 

\ht) = exp (-(2L_2 - (K/2)L^)ti) \ht-t^) • (53) 



But, as we argued earher, \ht) is independent of t, and therefore 

(2L_2- (/t/2)L2_,)|/i) =0. (54) 

This means that the descendant states L_2\h) and L'^_^\h) are linearly dependent. We say 
that the Virasoro representation corresponding to \h) has a null state at level 2. From 
this follow an number of important consequences. Acting on with Li and L2, and 
using the Virasoro algebra PHj) and the fact that Li\h) = L2\h) = {] while LQ\h) = h\h), 
leads to 

h = /i2,i = ; (55) 

c = <^^-«"°-^). (56) 

These are the fundamental relations between the parameter k, of SLE and the data of 
CFT. The conventional notation h2^i comes from the Kac formula in CFT which we do 
not need here. In fact this is appropriate to the case k < 4: for k > 4 it corresponds 
to hi^2- (To further confuse the matter, some authors reverse the labels.) Note that the 
boundary exponent h parametrises the failure of locality in (j2SI)- From CFT we may 
also deduce that, with respect to the path integral measure, the probability that a curve 
connects small intervals of size e about points ri, r2 on the real axis behaves like 

C\ 2/12,1 
I ' I . (57) 

Such a result, elementary in CFT, is difficult to obtain directly from SLE in which the 
curves are conditioned to begin and end at given points. 

Note that the central charge c vanishes when either locality {k = 6) or restriction (k = |) 
hold. These cases correspond to the continuum limit of percolation and self-avoiding 
walks respectively, corresponding to formal limits Q — ^ 1 in the Potts model and n — > 
in the 0{n) model for which the unconditioned partition function is trivial. 



5.4 Differential Equations 

In this section we discuss how the linear second order differential equations for various 
observables which arise from the stochastic aspect of SLE follow equivalently from the 
null state condition in CFT. In this context they are known as the BPZ equations As 
an example consider Schramm's formula (j30|l for the probability P that a point ( lies to 
the right of 7, or equivalently the expectation value of the indicator function 0{() which 
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is 1 if this is satisfied and zero otherwise. In SLE, this expectation value is with respect 
to the measure on curves which connect the point to infinity. In CFT, as explained 
above, we can only consider curves which intersect some e-neighbourhood on the real axis. 
Therefore P should be written as a ratio of expectation values with respect to the CFT 
measure: 

P(C;ao)^...n <^"'°°'°K'^"^;-» . (58) 

(02,1(00)02,1(^^2)) 

We can derive differential equations for the correlators in the numerator and denomina- 
tor by inserting into each of them a factor (l/(27ri) a{z)T{z)dz+ c.c, where a{z) = 
2/{z — ao), and F is a small semicircle surrounding ao- This is equivalent to making the 
infinitesimal transformation z z + 2e/{z — ao). As before, the c.c. term is equivalent to 
extending the contour in the first term to a full circle. The effect of this insertion may be 
evaluated in two ways: by shrinking the contour onto ao and using the OPE between T 
and 02,1 we get 

2L_202,i(ao) = {K/2)L\<P2,iiao) = d^Mao) , (59) 
while wrapping it about ( (in a clockwise sense) we get 

- (2/(C - ao))d^O - (2/(C - ao))d^O . (60) 

The effect on 02,1(^2) vanishes in the limit r2 00. As a result we can ignore the variation 
of the denominator in this case. Equating and inside the correlation function in 
the numerator then leads to the differential equation (j^^ for P found in Sec. 14.11 



5.4.1 Calagero-Sutherland model 

While many of the results of SLE may be re-derived in CFT with less rigour but perhaps 
greater simplicity, the latter contains much more information which is not immediately 
apparent from the SLE perspective. For example, one may consider correlation functions 
(0i,2('^i)0i,2('^2) • • •02,i('^Af) • • •) of multiple boundary condition changing operators with 
other operators either in the bulk or on the boundary. Evaluating the effect of an inser- 
tion (l/27ri) Jj.T{z)dz/{z — Vj) where F surrounds rj leads to a second order differential 
equation satisfied by the correlation function for each j. 

This property is very powerful in the radial version. Consider the correlation function 

C$(^i, . . . , Om) = (02,i(^i) . . . M^nMO)) (61) 

of 02.1 operators on the boundary of the unit disc with a single bulk operator $ at the 
origin. Consider the effect of inserting (l/27ri) aj{z)T{z)dz where (cf. (EHl)) 

, , z + e^'^^ , , 

and F surrounds the origin. Once again, this may be evaluated in two ways: by taking 
F up to the boundary, with exception of small semicircles around the points e*^*, we get 
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GjCq,, where Gj is the second order differential operator 

r =--—^^^ — -S"( ~ ^ ^ 

2dey 6 +12 2 36, 2sin\eu-e,)/2 

The first three terms come from evaluating the contour integral near e*^^, where Uj acts 
like 2L_2 — \Lq — j^, (the term ^ comes from the curvature of the boundary,) and the 
term with k j from the contour near e*^'-', where it acts like aj(e*^*)L_i + Rea^(e*^'=)Lo. 

On the other hand, shrinking the contour down on the origin we see that aj{z) = z+0{z^), 
so that on $(0) it has the effect of Lq + + . . ., where the omitted terms involve the L„ 
and Ln with n > 0. Assuming that $ is primary, these other terms vanish, leaving simply 
(Lq + ivo)$ = Equating the two evaluations we find the differential equation 

Gj C$ = a;$ C$ . (64) 




In general there is an {N — l)-dimensional space of independent differential operators 
Gj with a common eigenfunction C$. (There is one fewer dimension because they all 
commute with the total angular momentum J2j{d/d6j).) For the case N = 2, setting 
9 = 02 — Oi, we recognise the differential operator in Sec. 14.3 .31 

In general these operators are not self-adjoint and their spectrum is difficult to analyse. 
However, if we form the equally- weighted linear combination G = J2f=i Gj, the terms with 
a single derivative may be written in the form J2k{9V/d6k){d/d6k) where is a potential 
function. In this case it is well known from the theory of the Fokker-Plank equation 
that G is related by a similarity transformation to a self-adjoint operator. In factjSEl 
if we form \^ iy\'^/'^G\'^ n\~'^^'^ where = I\j<k{^^^^ ~ ^^^) is the 'free-fermion' wave 
function on the circle, the result is, up to calculable constants the well-known A^-particle 
Calogero-Sutherland hamiltonian 

H iB^- If I /^(/^-2) ^ 1 .... 

-"2^5^ + g sin^(^, - e,)/2 ' 

with f3 = 8/k,. It follows that the scaling dimensions of bulk operators like $ are simply 
related to eigenvalues A^v of by 

x$ = {k/N)An{8/k) - {A/kN)EI + |/i2,i + j^c , (66) 

where = ^N{N'^ — 1). Similarly C$ is proportional to the corresponding eigen- 
function. In fact the ground state (with conventional boundary conditions) turns out 
to correspond to the bulk A^-leg operator discussed in Sec. 12.4.21 The corresponding 
correlator is |\l/ArP/'^. 
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6 Related ideas 

6.0.2 Multiple SLEs 

We pointed out earlier that the boundary operators 02, i correspond to the continuum 
limits of lattice curves which hit the boundary at a given point. For a single curve, 
these are described by SLE, and we have shown in that case how the resulting differential 
equations also appear in CFT. Using the iV-particle generalisation of the CFT results 
of the previous section, we may now 'reverse engineer' the problem and conjecture the 
generalisation of SLE to curves. 

The expectation value of some observable O given that curves, starting at the origin, 
hit the boundary at {9i, . . . ,9]y) is 

Po(^i,...,^^) = — — , (67) 

J^l{t/i, . . . ,t/N) 

where Fq = (C02,i(e'^i) . . . 02,i(e^^^)<l>7v(O)). This satisfies the BPZ equation 

G, Fa = ((5,O02,i(e*'^) . . . 02,i(e^'")$ ^(0)) , (68) 

where 5jO is the variation in O under aj. If we now write Fq = Fi ■ Pq and use the 
fact that GjFi = we find a relatively simple differential equation for Pq, since the 

non- derivative terms in Gj cancel. There is also a complication since the second derivative 
gives a cross term proportional to {dg.Fi){d0^Po)- However, this may be evaluated from 
the explicit form F^ = |\l/7vp^'^. The result is 




Po = S,Po, (69) 



where the right hand side comes from the variation in O. 

The left hand side may be recognised as the generator (the adjoint of the Fokker-Planck 
operator) for the stochastic process 

dOj = ^dBt + cot {{9 J- 9k) /2)dt; (70) 

d9k = cot{{9k-9j)/2)dt, (71) 

where pk = 2. [For general values of the parameters pk this process is known as (radial) 
SLE(fi;,p), although this is more usually considered in the chordal version. It has been 
argued [37j that this applies to the level lines of a free gaussian field with piecewise constant 
Dirichlet boundary conditions: the parameters pk are related to the size of the disconti- 
nuities at the points e*^*. SLE(k, p) has also been used to give examples of restriction 
measures on curves which are not reflection symmetric! 
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We see that e*^-' undergoes Brownian motion but is also repelled by the other particles 
at e^^'' {k 7^ j): these particles are themselves repelled deterministically from e*^^. The 
infinitesimal transformation aj corresponds to the radial Loewner equation 

dgj,t _ gj,t + e*^^-* , . 

The conjectured interpretation of this is as follows: we have non-intersecting curves 
connecting the boundary points e*^*'° to the origin. The evolution of the jth curve in the 
presence of the others is given by the radial Loewner equation with, however, the driving 
term not being simple Brownian motion but instead the more complicated process ()70I7H) . 

However, from the CFT point of view we may equally well consider the linear combination 
J2j Gj- The Loewner equation is now 

^ gt + e*^^'* 

where 

dOj = ^dBi + 2 ^ cot {{Oj - Ok)/ 2) dt . (74) 

This is known in the theory of random matrices as Dyson's Brownian motion. It describes 
the statistics of the eigenvalues of unitary matrices. The conjectured interpretation is now 
in terms of random curves which are all growing in each other's mutual presence at 
the same mean rate (measured in Loewner time). From the point of view of SLE, it is 
by no means obvious that the measure on curves generated by process ()70l71l72j) is 
the same as that given by ()73I74|) . However CFT suggests that, for curves which are the 
continuum limit of suitable lattice models, this is indeed the case. 



6.1 Other variants of SLE 

So far we have discussed only chordal SLE, which describes curves connecting distinct 
points on the boundary of a simple connected domain, and radial SLE, in which the curve 
connects a boundary point to an interior point. Another simple variant is dipolar SLE|39j. 
in which the curve is constrained to start at boundary point and to end on some finite 
segment of the boundary not containing the point. The canonical domain is an infinitely 
long strip, with the curve starting a point on one edge and ending on the other edge. This 
set-up allows the computation of several interesting physical quantities. 

The study of SLE in multiply-connected domains is very interesting. Their conformal 
classes are characterised by a set of moduli, which change as the curve grows. R. Friedrich 
and co-workers jlO] have argued that SLE in such a domain is characterised by diffusion 
in moduli space as well as diffusion on the boundary. 
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It is possible to rewrite the differential equations which arise from null state conditions in 
extended CFTs (for example super-conformal CFTsjUj and WZWN models in terms 
of the generators of stochastic conformal mappings which generalise that of Loewner. 
However, a physical interpretation in terms of the continuum limit of lattice curves appears 
so far to be missing. 

6.2 Other growth models 

SLE is in fact just one very special, solvable, example of an approach to growth processes 
in two dimensions using conformal mappings which has been around for a number of 
years. For a recent review see The prototypical problem of this type is diffusion- 
limited aggregation (DLA). In this model of cluster formation, particles of finite radius 
diffuse in, one by one, from infinity until they hit the existing cluster, where they stick. 
The probability of sticking at a given point is proportional to the local electric field, if 
we imagine the cluster as being charged. The resultant highly branched structures are 
very similar to those observed in smoke particles, and in viscous fingering experiments 
where one fluid is forced into another in which it is immiscible. Hastings and LevitoV(44jj 
proposed an approach to this problem using conformal mappings. At each time t, the 
boundary of the cluster is described by the conformal mapping ft{z) which takes it to the 
unit disc. The cluster is grown by adding a small semicircular piece to the boundary. The 
way this changes ft is well known according to a theorem of Hadamard. The difficulty is 
that the probability of adding this piece at a given point depends on the local electric field 
which itself depends on The equation of motion for ft is therefore more complicated 
than in SLE. Moreover it may be shown that almost all initially smooth boundary curves 
evolve towards a finite-time singularity: this is thought to be responsible for branching, 
but just at this point the equations must be regularised to reflect the finite size of the 
particles (or, in viscous fingering, the effects of finite surface tension.) 

It is also possible to generate branching structures by making the driving term at in 
Loewner's equation discontinuous, for example taking it to be a Levy process. Unfortu- 
nately this does not appear to describe a physically interesting model. 

Finally, Hastings |45j has proposed two related growth models which each lead, in the 
continuum limit, to SLE. These are very similar to DLA, except that growth is only 
allowed at the tip. The first, called the arbitrary Laplacian random walk, takes place on 
the lattice. The tip moves to one of the neighbouring unoccupied sites r with relative 
probability E{r)^, where E{r) is the lattice electric field, that is the potential difference 
between the tip and r, and 77 is a parameter. The second growth model takes place in the 
continuum via iterated conformal mappings, in which pieces of length ii are added to the 
tip, but shifted to the left or right relative to the previous growth direction by a random 
amount ±£2- This model depends on the ratio £2/^1, and leads, in the continuum limit, 
to SLEk with K = £2/4£i- For the lattice model there is no universal relation between k 
and 1], except for t] = 1, which is the same as the loop-erased random walk (Sec. 12. 2|) and 
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converges to SLE2. 
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